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PACS numbers: 13.30.-a 13.30.Eg 13.40.-f 13.60.-r 13.60.Le 14.20.Gk 25.20.Lj 



I. INTRODUCTION 

We report on measurements of the photoproduction 
from the proton of two ground state hyperons, namely 
the reactions 7 + p — » K + + A and 7 + p — > K + + T°. In- 
termediate baryonic states in these reactions can be the 
N* resonances in the case of A production, and N* or A 
resonances in the case of T° production. In either case 
one expects strange meson exchange in the t channel and 
hyperon exchange in the u channel. This is illustrated 
in Fig.^ To unravel the production mechanism in these 
reactions, highly detailed measurements of as many ob- 
servables as possible are needed. 

In this paper we present results for the differential cross 
sections, da/ d cos obtained with the CLAS sys- 
tem in Hall B at Jefferson Lab. Following our previous 
publication, Ref. pj, these results are based on additional 
data accumulated by CLAS and use a different analysis 
technique. In another forthcoming paper we will present 
results for the beam-recoil double polarization observ- 
ables, C' x and C z , for the same reactions obtained from 
the same data set. 

The main motivation for this work was to provide 
data to investigate the spectrum of non-strange (N* and 
A) baryon resonances above the strangeness-production 
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threshold at W = -y/s = 1.6 GeV. Between this threshold 
and the upper limit of our data set, at W = 2.53 GeV, 
many baryon resonances are predicted by quark mod- 
els 2] , but relatively few are clearly established 3] . These 
resonances are broad and overlapping, making partial 
wave analysis challenging, but it is also possible that 
some dynamical aspect of hadronic structure may act to 
restrict the quark models' spectrum of states to some- 
thing closer to what has already been established Q. 
This is the so-called "missing resonance" problem. While 
the branching fractions of most high-mass resonances to 
K Y final states are expected to be small (cross sections 
~ 1 fib) compared to three-body modes such as 7T7t N 
(~ 100 fib), the study of these decays do have advan- 
tages. First, two-body final states are often easier to 
analyze than three-body final states. Second, couplings 
of nucleon resonances to KY final states will differ from 
coupling to ttN, rjN, or nnN final states 0. Thus, one 
can hope that this alternate light cast on the baryon reso- 
nance spectrum may emphasize resonances not otherwise 
revealed. Some "missing" resonances may only be "hid- 
den" when sought in more well-studied reaction channels. 

The A and S° hyperons have isospin and 1, re- 
spectively, and so intermediate baryonic states leading 
to the production of A's can only have isospin 1/2 (TV* 
only), whereas for the TP's, intermediate states with both 
isospin 1/2 and 3/2 (N* or A) can contribute. Thus, si- 
multaneous study of these reactions provides a kind of 
isospin selectivity of the sort used in comparing 77 and 7r 
photoproduction reactions. To date, however, the PDG 
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FIG. 1: Representative tree-level diagrams illustrating s— 
(top), t— (middle), and u— channel (bottom) exchanges. 
Born terms (left column), baryon resonance excitations (top 
right), and other exchanges (right middle, bottom) lead to 
production of K Y . Models differ in their electromagnetic 
transition moments (/x's),the strong couplings (gMBs's), and 
form factors, as well as the effects of channel couplings. 



compilation ,3] gives poorly-known KA couplings for only 
five well-established resonances, and no ifS couplings for 
any resonances. The most widely-available model cal- 
culation of the KA photoproduction, the Kaon-MAID 
code [f|, includes a mere three well-established N* states: 
the 5n (1650), the Pn (1710), and the P 13 ( 1720). Thus, it 
is timely and interesting to have additional good-quality 
photoproduction data of these channels to see what ad- 
ditional resonance formation and decay information can 
be obtained. 

Section II of this paper discusses briefly the reaction 
models that will be compared with the present data. Sec- 
tion III discusses the experimental setup of the CLAS sys- 
tem for this experiment. The steps taken to obtain the 
cross sections from the raw data are discussed in Section 
IV. Section V presents the results for the measured angu- 
lar distributions and VP-dependence of the cross sections. 
In Section VI we discuss the results in light of previ- 
ous measurements, and in relation to several previously- 
published reaction models. We also show how the data 
can be parameterized in terms of t-channel scaling us- 
ing a simple Regge-based picture, and in terms of simple 
Legendre polynomials. In Section VII we recapitulate the 
main results. 



II. THEORETICAL MODELS 

The results in this experiment will be compared to 
model calculations that fall into two classes: tree- level 
Effective Lagrangian models and Reggeized meson ex- 
change models. Effective Lagrangian models evaluate 
tree-level Feynman diagrams as in Fig. ^ including reso- 
nant and non-resonant exchanges of baryons and mesons. 
A complete description of the physics processes will re- 
quire taking into account all possible channels which 
could couple to the one being measured, but the advan- 
tages of the tree-level approach are to limit complex- 
ity and to identify the dominant trends. In the one- 
channel tree-level approach, some tens of parameters (in 
particular, the couplings of the non-strange baryon res- 
onances to the hyperon-kaon systems) must be fixed by 
fitting to data, since they are poorly known from other 
sources. An alternative approach is to use no baryon 
resonance terms and instead model the cross sections in 
a Reggeized meson exchange picture. While this is not 
expected to reproduce the results in detail, it will show 
where the high-energy phenomenology of i-channel dom- 
inance blends into the nucleon resonance region picture. 

For K + A production, the model of Mart and 
Bennhold |(| has four baryon resonance contributions. 
Near threshold, the steep rise of the cross section is ac- 
counted for with the N* states S , n(1650), the P u (1710), 
and -Pi3(1720). To explain the broad cross-section bump 
in the mass range above these resonances, they intro- 
duced the Z?i3(1895) resonance that was predicted in the 
relativized quark models of Capstick and Roberts [2j and 
Loring, Metsch, and Petry to have especially strong 
coupling to the K + A channel. In addition, the higher 
mass region has contributions, in this model, from the 
exchange of vector K*(892) and pseudovector ifi(1270) 
mesons. The hadronic form factors, cutoff masses, and 
the prescription for enforcing gauge invariance were ele- 
ments of the model for which specific choices were made. 
The content of this model is embedded in the Kaon- 
MAID code |5j which was used for the comparisons in 
this paper. This model was tuned to results from the 
experiment at Bonn/SAPHIR and offers a fair de- 
scription of those results. 

On the other hand, analysis by Saghai et al. |9( using 
the same data set showed that, by tuning the background 
processes involved, the need for the extra resonance was 
removed. Janssen et al. |ToL ITlj showed that the same 
data set was not complete enough to make firm state- 
ments since models with and without the presence of a 
hypothesized iV*(1895)Di3 resulted in equally good fits 
to data. A subsequent analysis [l2|, which also fitted 
calculations to photon beam asymmetry measurements 
from SPring-8 Il3| and electroproduction data measured 
at Jefferson Lab [l4( , indicated weak evidence for one or 
more of S\\, P\\, P13, or Z?i 3 (1895), with the P\\ solution 
giving the best fit. The conclusion was that a more com- 
prehensive data set would be required to make further 
progress. 
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More elaborate model calculations have been under- 
taken in which channel coupling is considered, in addition 
to the tree-level approaches mentioned above. Penner 
and Mosel found fair agreement for the K + A data 
without invoking a new D13 structure. Chiang et al. [l6l | 
showed that coupled channel effects are significant at the 
20% level in the total cross sections when including pio- 
nic final states. Shklyar, Lenske, and Mosel [Tjj used a 
unitary coupled-channel effective Lagrangian model ap- 
plied to 7r and 7 -induced reactions to find dominant 
resonant contributions from Sii(1650), Pi3(1720), and 
Pi 3 (1895) states, but not from P n (1710) or £> 13 (1895). 
This conclusion was true despite the discrepancies be- 
tween previous data from CLAS and SAPHIR (Isj . 
Recently, Sarantsev et al. did a phenomenological 
multi-channel fit for KA, KT,, as well as tt and 77 photo- 
production data. They found fairly strong evidence for 
a P n at 1840 MeV and two D 13 states at 1870 and 2170 
MeV. Even better quality KY data such as we are pre- 
senting here are needed to solidify these conclusions. We 
will not compare the present results to those models in 
this paper, however. 

While it is to be expected that s-channel resonance 
structure is a significant component of the K + A and 
K + Ti° reaction mechanisms, it is instructive to compare 
to a model that has no such content at all. The model 
of Guidal, Laget, and Vanderhaeghen is such a 

model, in which the exchanges are restricted to two lin- 
ear Regge trajectories corresponding to the vector K* 
and the pseudovector K\. The model was fit to higher- 
energy photoproduction data where there is little doubt 
of the dominance of these exchanges. In this paper, we 
extend that model into the resonance region in order to 
make a critical comparison. 



III. EXPERIMENTAL SETUP 

Differential cross section data were obtained with the 
CLAS system in Hall B at the Thomas Jefferson Na- 
tional Accelerator Facility. Electron beam energies of 2.4 
and 3.1 GeV contributed to the data set, each of typ- 
ically 10 nA current. Real photons were produced via 
bremsstrahlung from a 1 x 10~ 4 radiation length gold 
radiator and "tagged" using the recoiling electrons ana- 
lyzed in a dipole magnet and scintillator hodoscopes (2^ • 
The energy tagging range was from 20% to 95% of the 
beam endpoint energy, and the integrated rate of tagged 
photons was typically 5 x 10 6 /sec. Using the tagger and 
the accelerator RF signal, photon timing at the physics 
target was defined with an rms precision of 180 psec. 
The useful energy range for this experiment was from 
the strangeness-production threshold at E~, — 0.911 GeV 
(W = 1.61 GeV) up to 2.95 GeV (W = 2.53 GeV). In 
this range, the tagger resolution was typically 5 MeV, set 
by the size of the hodoscope elements, but the data were 
analyzed in bins of 25 MeV photon energy to be commen- 
surate with any energy-dependent structure expected in 



the hadronic cross sections. The centroids of these bins 
were adjusted in the analysis by between —6 and +5 MeV 
to compensate for mechanical sag of the hodoscope ar- 
ray measured by kinematically fitting p(-y,pTT + TT~) data; 
hence our final results are given in unequal energy steps. 

The physics target consisted of a 17.9 cm long liquid 
hydrogen cell of diameter 4.0 cm. Temperature and pres- 
sure were monitored continuously to determine the den- 
sity to 0.3% precision. The target cell was surrounded 
by a set of six 3 mm thick scintillators to help define 
the starting time for particle tracks leaving the target, 
though actually the timing given by the photon tagger 
was used to define the event times. 

The CLAS system, described in detail elsewhere [23| . 
consisted of a toroidal magnetic field, with drift chamber 
tracking of charged particles. The overall geometry was 
six-fold symmetric viewed along the beam line. Particles 
could be tracked from 8° to 140° in laboratory polar an- 
gle, and over about 80% of 2ir in the azimuthal direction. 
Outside the magnetic field region a set of 288 scintillators 
was used for triggering and for later particle identifica- 
tion using the timc-of-flight technique. The momentum 
resolution of the system was ~ 0.5%, with variations due 
to multiple scattering and tracking resolution considera- 
tions. The low-momentum cut-off was set in the analysis 
at 200 MeV/c. Other components of the CLAS system, 
such as the electromagnetic calorimeter and the Cerenkov 
counters, were not used for these measurements. 

The event trigger required an electron signal from the 
photon tagger, and at least one charged-track coincidence 
between the time-of- flight 'Start' counters near the tar- 
get and the time-of-flight 'Stop' counters surrounding the 
drift chambers. The photon tagger signal consisted of 
the OR of coincidences among hits in a two-plane ho- 
doscope, which had 61 timing scintillators in coincidence 
with their matching energy-defining scintillators. The 
charged-track trigger in CLAS was a coincidence of six 
OR'd start counter elements and the OR of the outer 
time-of-flight scintillators. Events were accumulated at 
the rate of ~ 2500 hadronic events per second, though 
only a sub-percent fraction of these events contained the 
kaons and hyperons of interest for the present analysis. 



IV. DATA ANALYSIS 

A. Data and event selection 

The data used in this experiment were obtained in late 
1999 as part of the CLAS "gle" data taking period. Since 
the electronic trigger was loose, data for several photo- 
production studies were contained in the data set. Off- 
line calibration was performed to align the timing spec- 
tra of the elements of the photon tagger, the six elements 
of the start counter, and the 288 elements of the time- 
of-flight (TOF) counters. Drift-time calibrations were 
made for the 18 drift chamber packages. Pulse height 
calibrations and timing-walk corrections were made for 
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the time-of-flight counters. The raw data were then pro- 
cessed to reconstruct tracks in the drift chambers and to 
associate them with hits in the time-of-flight counters. 

B. Particle Identification 

Kaon, proton, and pion tracks were separated using 
momentum and time-of-fiight measurements. The mo- 
mentum, p, of each track was measured directly via track 
reconstruction through the CLAS magnetic field; this 
measurement also gave the path length, d, from the reac- 
tion vertex to the time-of-flight counter hit by the track. 
The starting time of the track was determined by pro- 
jecting the tagger signal time, synchronized with the ac- 
celerator RF timing, to the reaction vertex inside the 
hydrogen target. The stopping time was determined by 
the element hit in the array of TOF scintillators. The 
difference, T, between these two times was the measured 
time of flight, which in CLAS could range between about 
4 and 100 nsec. From T the speed, j3, could be obtained 
as (3 = d/(cT). The mass, m„ was then computed ac- 
cording to m, = \J\ — (3 2 x pT/d. In CLAS, the dom- 
inant mass uncertainty in this situation came from the 
time-of-flight resolution, ST. ST was independent of par- 
ticle momentum, so particle selection based on time of 
flight was largely independent of momentum as well. For 
kaon identification we used the time-of-flight difference 
technique, where the measured time, T, of the track was 
compared to the expected time, Th, for a hadron of mass 
rrih and momentum p. For a hypothesized value of mh 
we can define Atof = T — Th and write 

y y (m x cy +p 2 j 

Figure [21 shows an example of such a time difference 
spectrum when we took nih to be the kaon mass. The 
candidate kaon tracks were selected using a ±1 nsec cut 
centered at zero. Pion and proton bands are well sepa- 
rated from the kaons up to 1 and 2 GeV, respectively. A 
crossing band due to a badly-calibrated detector element 
is shown for illustration; such tracks were later rejected 
by removing the detector element and/or by the kine- 
matic cuts and fits applied later. Above 1 GeV some pi- 
ons leak into the set of candidate kaon tracks. These were 
rejected by subsequent event reconstruction cuts and by 
background rejection fitting. Protons were identified us- 
ing a similar Atof correlation but with looser cuts due to 
the straggling effects which broadened the distribution. 

Photons matching the hadronic tracks in CLAS were 
selected using the time difference between the hadronic 
track projected back to the event vertex and the photon 
tagger time projected forward to the event vertex. Fig- 
ure |3 shows such a spectrum, which illustrates the pres- 
ence of random coincidences between the photons and 
the hadronic tracks. The 2 nsec RF time structure of the 
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FIG. 2: (Color online) Time-of-flight difference spectrum for 
a sample of tracks, assuming the mass of the particle is that 
of a kaon. White lines indicate the cut limits for selecting 
kaons in a time window of ±1.0 nsec. Note the logarithmic 
scale on the intensity axis. 
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FIG. 3: Time difference between photon tagger time and 
target start counter time showing the peak at zero of good 
matches between the photons and the hadrons at the event 
vertex in CLAS. Coincidences due to hadrons mismatched to 
random photons in the tagger show the 2 nsec bunch structure 
of CEBAF. 

accelerator is clearly seen. A ±1.0 nsec cut was used to 
reject out-of-time combinations. In-time accidentals un- 
der the central peak were treated as potentially-correct 
photons, and such particle-photon combinations were re- 
tained in the analysis. Since ambiguous photons were 
generally widely separated in energy, the (7, K + ) miss- 
ing mass for incorrect combinations fell into the broad 
background under the hyperons, and were then rejected 
at the peak-fitting stage of the analysis discussed below. 

In this analysis we demanded detection of positive 
kaons and protons. Negative pions from A decay or pho- 
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tons from S° decay were not required. Fiducial cuts were 
applied in track angle and momentum to restrict events 
to the well-described portions of the detector. This in- 
cluded removal of 9 out of 288 time-of-flight elements 
due to poor timing properties. Corrections were applied 
for the mean energy losses of kaons and protons as they 
passed through the production target, target walls, beam 
pipe, and air. The nominal CLAS momentum reconstruc- 
tion algorithms were found to provide sufficient hyperon 
mass resolution (see below) that no higher-order momen- 
tum corrections were applied. 

A missing mass cut was applied to p(~f, K + p)ir~ (7) 
to select events consistent with a missing pion and (for 
the £ ) a missing photon. The losses incurred by this cut 
due to multiple scattering effects on the part of the kaons 
and protons were studied in the real data and in Monte 
Carlo. The estimated residual uncertainty due to the cut 
and its compensation via the acceptance calculation was 
l%-2%. 
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FIG. 4: Hyperon spectrum via missing mass using the photon 
and detected kaon, integrated over all kaon angles and photon 
energies using a 3.1 GeV endpoint energy. 



C. Yield of hyperons 

The extraction of kaon yields in each bin of photon 
energy and kaon angle depended on fits to the missing 
mass spectrum given by p("f, K + )Y . When integrating 
over all of our 3.1 GeV data, for all energies and angles, 
the resulting missing mass spectrum is shown in Fig. H 
This figure illustrates that the overall missing mass res- 
olution of the system was a — 8.9 MeV for the A and 
a = 8.2 MeV for the E°. The overall resolution averaged 
6.3 MeV in the 2.4 GeV data set, where all the average 
momenta were lower. However, the width of the peaks 
and the extent of the background to be removed from un- 
der the peaks via fitting varied substantially across the 
measured range of energy and angle, so a careful fitting 
procedure was needed to obtain well-controlled hyperon 
yields. 

The main source of background in the hyperon mass 
spectra was due to events where the kaons were actually 
mis-identified pions. The yields of A and E° hyperons 
were obtained using lineshape fits to missing-mass spec- 
tra in each of over 1,450 kinematic bins of photon energy 
and kaon angle. The data were binned in 25 MeV steps 
in E~f and in 18 bins of kaon center-of-mass (cm.) an- 
gle, cos(#^-™), centered in steps of 0.1 between —0.8 and 
+0.9. 

Typical hyperon yield fits of p(j, K + )Y for the middle 
of the photon energy range are shown in Fig. |3J The fits 
were performed in two passes. In the first pass, events 
for all kaon center-of-mass angles were summed together. 
These first fits served to determine and fix the centroids 
and widths of the Gaussian peaks for the two hyperons. 
These were 7 to 9 parameter fits, depending on the back- 
ground model employed. A log-likelihood fitting algo- 
rithm was used. The background was modeled as poly- 
nomials of order up to 2 (quadratic). In the second pass, 
fits were made with 3 to 5 parameters for the yields in 



each kaon angle bin, allowing only the integrated counts 
of the peaks to vary in addition to the background pa- 
rameters. The two-pass method was used to stabilize fits 
of low-yield bins at low photon energy and backward kaon 
angle. Acceptable fits all had \ 2 P er degree of freedom 
of less than 2.0. 

Background parameterizations that were a simple con- 
stant or a sloped line were sufficient to yield good fits over 
most of the kinematic range. At more forward kaon an- 
gles the effect of background due to mis-identified pions 
increased and the quadratic fits generally gave the best 
results. Above E 1 — 2.3 GeV the momentum resolution 
of CLAS broadened to the degree that the forward-angle 
quadratic fits became less stable, so the linear fits were 
preferred. This led to an extra estimated systematic un- 
certainty of 10% on both the forward-angle differential 
and total cross sections above this energy. In some low- 
yield back-angle bins, where no good fits were obtained, 
side-band subtraction was used to determine the yield. 

The final cross sections were based on the following 
numbers of fully reconstructed events: from the 2.4 and 
the 3.1 GeV data sets we had 236,260 and 325,792 K+A 
events, respectively, and 169,796 and 269,216 A+£ 
events, respectively. 



D. Acceptance calculation 

The acceptance and efficiency were modeled us- 
ing a CLAS-standard GEANT-based simulation code 
("GSIM"). Events were initially modeled using a phase 
space distribution for 7 + p -> K+ + Y. The GSIM 
code simulated the events in the CLAS detector at the 
level of ADC and TDC hits in the scintillators and drift 
time information in the tracking chambers. The events 
were further fine-tuned such that the time distributions 
in the TDC's accurately matched the actual data us- 
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FIG. 5: Sample missing mass fits used for the determination of hyperon yields at E-y — 1.825 GeV and for three representative 
K + angles: cos(6^™') = —0.7, +0.1, +0.8. The fits to different orders of polynomial background are nearly indistinguishable: 
solid lines for quadratic and dashed lines for linear. 



ing another well-tested CLAS software package ("GPP"). 
These simulated events were then processed through the 
same analysis codes as the real data, and thus the ac- 
ceptance was computed in each kinematic bin. Dead 
regions of the drift chambers were removed both from 
the real data and from the simulated data during track 
reconstruction ("A1C"). Detector efficiency was simulta- 
neously accounted for through the simulation: sources 
of inefficiency included track reconstruction failures and 
time-of- flight paddle removals. The only particle back- 
ground in this physics Monte Carlo was due to particle 
decays, especially the kaons, and multiple scattering ef- 
fects. Thus, we relied on the yield extractions discussed 
earlier to remove background due to mis-identified pions 
or protons. 

The effect of using a phase-space event generator to 
compute the acceptance, rjp.s., was studied by using the 
fits to the angular distributions presented in Section IVAl 
to regenerate the acceptance, r\Bata-, with an improved 
representation of the reactions. Since these cross sec- 
tions vary quite slowly with angle, and since the kine- 
matic bins were each small on the scale of these varia- 
tions, no large effects were to be expected. We found 
agreement between the two acceptance models at the 
level of 0.25% rms over essentially the whole of the kine- 
matic space, consistent with the statistical variations of 
the simulations. The exception was in the forward- most 
angle bin (0.85 < cos(<?^?) < 0.95) for both hyper- 
ons. There, because of the extrapolation of the anal- 
ysis into CLAS's forward acceptance hole, the ratio of 
acceptances rj£, ata /rip.s. dropped from 1.0 to 0.85 over 
the range £ 7 - 1.75 GeV to E 1 ~ 2.90 GeV. Theoreti- 
cal models of the behavior of the cross section in the very 
forward direction differ strongly, as shown later, so it was 
not known whether a "flat" or a "forward-peaked" or a 
"forward-dipped" acceptance model was more accurate. 
Thus, the forward-most angle results at cos(6'^ 1 " ) = +.9 
have an additional systematic uncertainty on the cross 
section which is, on the average, ±8%. 



K + A acceptance, 3. 1 GeV dataset 




cos(0 KCM ) 

FIG. 6: (Color online) Computed CLAS acceptance for the 
E-y — 3.1 GeV data set for the 7 + p — > K + + A reaction. The 
scale on the right gives the value of the acceptance for each 
kinematic bin. 



A sample of the acceptances computed for CLAS for 
these reactions is shown in Fig. [S] It was largest at mid 
to forward kaon angles and at higher photon energies. 
The maximum acceptance was about 22% for K + A and 
23% for K + YP . A lower cut-off was applied, such that 
the smallest allowed acceptance in the experiment was 
0.5%. For each hyperon, 10 million events were gener- 
ated at each beam endpoint energy. Non-uniformities in 
the distribution arise from the effects of detector element 
removals and track reconstruction efficiencies. Since the 
kinematics of the two hyperon reactions are very similar, 
the acceptance function for the E° looked very similar, 
apart from the higher production threshold. 
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E. Photon flux 

The number of photons striking the target was com- 
puted from the measured rate of electrons detected at 
the tagger hodoscope. TDC spectra of the tagger ele- 
ments recorded the hits of electrons in a 150 nsec time 
range around each event. This flux was scaled and in- 
tegrated in ten-second intervals. After statistical correc- 
tions for multiple hits and electronic live time, the flux 
was summed over whole runs. The fine granularity of 
the tagging system was grouped into bins of 25 MeV in 
photon energy. 

Photon losses in the beam line due to tagger accep- 
tance, beam collimation, and thin windows were deter- 
mined using a separate total-absorption counter down- 
stream of CLAS. This low-rate lead glass detector was 
periodically put in the beam line to monitor the tagging 
efficiency. For the 2.4 GeV data set the average efficiency 
for tagging photons was 78%, and the stability of this ef- 
ficiency, which was measured periodically throughout the 
data taking period, was ±0.5%. 

By taking data at 2.4 and 3.1 GeV endpoint energies 
it was possible to test the flux normalization of many 
elements of the tagging system, as discussed in the next 
section. At energies above E 1 — 2.325 GeV the two 
data sets no longer overlapped, however, and defective 
electronics in a few channels of the tagger led to a gap 
in our final spectra. Bins at E 1 — 2.375 and 2.400 GeV 
were removed because of this. 



F. Systematic uncertainties 

The 2.4 and 3.1 GeV photon beam endpoint data sets 
were compared to investigate variations in the photon 
tagger efficiency. The photon-normalized yield of par- 
ticles at any given energy had to be independent of 
bremsstrahlung endpoint energy, so consistency of this 
quantity tested stability of the electronics. Localized re- 
gions of tagger inefficiency "moved" in photon energy 
when the endpoint energy changed. We took the higher 
normalized yield between the two data sets as the correct 
one. Localized regions of high inefficiency were found in 
the 3.1 GeV data set at 1.1, 1.4, and 1.8 GeV; in those 
regions we made corrections of up to 50% in one data set 
to compensate for tagger efficiency losses in the other. 
Much smaller corrections (~ 3%) were made at other en- 
ergies. The absolute uncertainty on these corrections was 
estimated to be ±3%. 

As a check on our results, the p(-f,ir + )n cross sec- 
tion was measured using the same analysis chain, as far 
as possible, as the p(j,K + )Y data. The same proce- 
dure was also used to generate the acceptance for the 
p(7, n + )n cross sections used to check the whole analysis 
process, except that the SAID code was used to generate 
the initial events. Figure \7\ shows the pion cross section 
measured in this analysis as a function of W for a mid- 
range cm. pion angle. The CLAS pion cross section 
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FIG. 7: CLAS differential cross section as a function of W 
for 7 + p — * 7r + + n. Shown for comparison are two versions 
of the SAID parameterization of world pion data. They are 
essentially indistinguishable. The vertical dotted line is the 
strangeness production threshold. 

was found to be in fair agreement with the SAID [24| pa- 
rameterization of the world's data between W = 1.6 and 
2.1 GeV, albeit lower by an average scale factor of 0.95. 
As a function of pion center of mass angle, the CLAS 
to SAID ratio was ~ 1.0 at back angles and ~ .92 at 
forward angles. Thus the pion results indicate a possi- 
ble systematic error in the acceptance calculation at the 
level of ±3%, apart from the average scale factor. The 
absolute accuracy of the pion cross sections, as parame- 
terized by SAID over the range of comparison we used, 
is similar to this. Therefore, we chose not to make a 
renormalization of our results to the average pion cross 
sections. The results presented in this paper are on an 
absolute scale. The kaon analysis was not identical to the 
pion analysis, since the kaon decay corrections are much 
larger in the former case, and since the final kaon anal- 
ysis included detection of the proton from the hyperon 
decays. Hence, it was difficult to translate the systematic 
trends in our pion results compared to SAID to the kaon 
results presented here. However, based on the compar- 
ison to the pion data analysis, we estimated the overall 
systematic uncertainty in our kaon cross section to be 
less than ±7%. This was the largest single contribution 
to our cross section systematic uncertainty. 

The analysis of this experiment was done twice, in 
largely independent ways. The first analysis |ll l25| com- 
puted cross sections based on detection of the kaon alone, 
or p(-f, K + )Y . Starting from the same data set, the sec- 
ond analysis [2|J detected the kaon and the proton for 
each event, or p(j, K + p)?t~ ('-/), where the tt~ (and the 
possible 7 from E° decay) were ignored. Particle identi- 
fication and acceptances were developed independently. 
For the results presented here, the first analysis was re- 
vised to take into account more advanced modeling of the 
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CLAS detector in the acceptance; both analyses used the 
standard CLAS GEANT package for computing accep- 
tances. The same flux normalization procedures were 
used. The first analysis used only data from the 2.4 GeV 
endpoint data set, while the second analysis also included 
data from 3.1 GeV endpoint. Consistency checks were 
then made between the two analyses. Results for the fi- 
nal cross sections from the two studies were in very good 
agreement across the full range of energies and angles 
where they overlapped. Isolated differences of ~ 5% in 
small ranges of angle were attributed to details of the 
acceptance modeling. By comparing the acceptances de- 
veloped over the course of the p(j, K + )Y studies, we es- 
timated that average systematic uncertainty across the 
kinematics of the experiment was ±2.7%, arising from 
variations in the implementation of the detector model 
and the track-reconstruction algorithms. 

On the energy axis, our results are precise to ±2 MeV. 
This systematic uncertainty arises from an energy bin 
centering correction that was applied to each data point 
due to the calibration of the photon tagger. In an 
independent study, kinematic fitting to the reaction 
p(7,p7r + 7r~) showed that the CLAS tagger and the pho- 
ton beam were mismatched by up to ±10 MeV due to 
mechanical effects in the structure of the tagger. The 
correction shifted the centroids of each energy bin by an 
amount estimated to be precise as stated above. The 
indirect effect that this centroid shift had on the accep- 
tance of CLAS was considered negligible, since the cross 
sections vary slowly in energy and the energy bins for the 
results are 25 MeV wide. 

The estimated systematic uncertainties discussed 
above were combined with contributions due to particle 
yield extraction (3.6%), photon attenuation in the beam 
line (0.2%), target density uncertainty (0.14%), and tar- 
get length uncertainty (0.28%). This led to an estimate 
of the global scale uncertainty of ±8%. Due to addi- 
tional systematic uncertainty about extrapolation of the 
data to zero degrees, the forward-most angle bin above 
E 1 = 1.75 GeV has an overall uncertainty of ±11%. 



V. RESULTS 

A. Angular Distributions 

Since the differential cross sections in this measure- 
ment are symmetric in the azimuthal angle 4>, we present 
the results in the partially integrated form 

da da 
d(cos 0^™-) = 71 dh ^ ' 

since this also puts the values on a convenient scale of 
order 1 /j.b. 

The angular distribution results for the reaction 7 + 
p — > K + +A are shown in Fig.[H] The results are shown as 
a function of cos(6^.™ ) for 79 bins in W. The step sizes in 



W were determined by the 25-MeV step size in the nom- 
inal photon energy, P 7 , at which the cross sections were 
extracted, together with a few-MeV correction for tag- 
ger re-calibration. There are 18 bins in cos(0S+'), each 
of width 0.1, centered from -0.80 to +0.90. The cross 
sections are the averages within each angle bin, with no 
bin centering. The results are the weighted means of the 
2.4 and 3.1 GeV beam energy data sets. The error bars 
are dominated by the statistical uncertainties of the hy- 
peron yield extraction fits, but also include the statistics 
from the Monte Carlo acceptances. The overall system- 
atic uncertainty, as discussed previously, is ±8%, except 
in the forward- most bin where above E 1 — 1.75 GeV it 
is ±11%. There are 1,377 data points in the K + A set. 

The curves in Fig. |H| arise from fits intended to capture 
the main features of the decay amplitudes contributing 
to the angular distributions. The form is 

where the Pi are the Legendre polynomials of order i, 
and the are the fit coefficients which represent the L — 
0, 1, 2,3,4= S, P, D, F, G -wave amplitudes for the decay 
distributions. The factor q/k is the phase space ratio of 
the reaction, where k and q are the center-of-mass frame 
momenta of the initial and final states, respectively. The 
value of this ratio ranges from zero at threshold up to .86 
at our highest energy. 

Qualitatively, the cross section is flat as a function of 
cos(0j^+') near threshold, as would be expected for S- 
wave behavior. As the energy rises to about 1.8 GeV 
the cross section develops a significant forward peaking 
consistent either with t-channel contributions or with s- 
channel interference effects between even and odd waves. 
As the energy rises further the cross section develops a 
tendency toward a slower rise in the extreme forward 
direction and also a rise in the backward direction. Above 
about 2.3 GeV the cross section is dominantly forward 
peaked, consistent with i-channel exchange dominance, 
though on a logarithmic scale (see discussion in Sec. I VI C|) 
the fall-off is not exponential all the way to back angles. 

The parameters of the fit may be used to gain some in- 
sight into the reaction mechanism, unraveling effects due 
to interference among partial waves. Figure El shows the 
coefficients from the fit using Eq. [21 The a 2 ; were taken 
to be purely real numbers. The range over which each 
parameter is plotted depended upon its significance, as 
estimated by the statistical F-test. Mostly, the higher 
partial waves are not significant near threshold, but our 
angular coverage is also less complete near threshold, dis- 
allowing higher-order fits. One may note a prominent 
bump in the P-wave amplitude between threshold and 
1.9 GeV, centered near 1.7 GeV. The D-wave amplitude 
turns on quite strongly near 1.9 GeV, and the P-wave 
amplitude has a broad dip centered at 2.05 GeV. In the 
K + A case, the G-wave was not significant at any energy. 

An alternative fitting procedure was performed that 
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FIG. 8: Differential cross sections for 7 + p — > K + + A. The number in each panel designates W (= 
results of the amplitude fits (Eq. [3J discussed in the text. 



y/s). The solid lines are 



decomposes the angular distribution magnitudes directly 
into Legendre coefficients, rather than amplitude-level 
partial wave Legendre coefficients. The fits were of the 
form 



da 



O 'tot 



dcos(8 c K n + l - 



4 

E 

i=l 



CiP,(cos0£?-) 



(4) 



and are shown in Fig. The total cross section, cr fot , 
was used as a parameter in order to obtain a proper es- 
timate of its uncertainty, which, due to parameter co- 
variances, is more difficult with the fits using Eq. [3J The 
coefficients are dimensionless ratios of the i th moments 
of the angular distribution to the total cross section. This 
fit procedure, to magnitudes rather than amplitudes of 
the distributions, is less useful in revealing interference 
effects. Nevertheless, some structure is visible. The C\ 
parameter shows a bump below 1.9 GeV which arises ei- 
ther from S-P or higher wave interference, and the C3 
parameter has a change in slope near 2.05 GeV. Over- 
all, the increasingly forward-peaked cross section with 
increasing energy forces all the CVs to rise with W. 



The differential cross sections for the A can be com- 
pared to the angular distributions for E° production 
shown in Fig. ^2 The bins in W are the same as before, 
allowing direct comparison of the panels in Figs. [SJand 
ITT1 Results for both hyperons were extracted together, 
using identical procedures discussed previously. There 
are 1,280 data points in the K + YP angular distributions. 

Besides the higher reaction threshold, the most signif- 
icant qualitative difference is that the E° cross section is 
not forward peaked in the energy range below 2 GeV. At 
W = 1.85 GeV, for example, the cross section peaks near 
cos(#^-™ 1 ') = 0.35, or 70° in the center-of-mass frame. 
This is consistent with a reaction mechanism for S° pro- 
duction that is less influenced by i-channel exchanges and 
is more s-channel resonance dominated than A produc- 
tion. The back-angle cross section is less prominent than 
for the A case in this energy range as well. Above the nu- 
cleon resonance region (above about 2.4 GeV), however, 
the two channels look quite similar, with characteristic 
i-channel forward peaking. 

The coefficients of the amplitude-level fit in Eq. [3] for 
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FIG. 9: (Color online) Amplitude fit to the differential cross sections for 7+p — > A" + + A. The coefficients are defined in Eq. [3] 
The solid vertical lines mark the well-known TV* resonances Sii(1650), Pn(1710), and Pi3(1720). The dotted line marks the 
E° threshold, and the dashed line marks the Di3(1895) position. 
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FIG. 10: (Color online) Fit to the magnitude of the differential cross sections for 7 + p — » K + + A. The coefficients are denned 
in Eq. 2] The solid vertical lines mark the well-known TV* resonances Sii(1650), Pn(1710), and Pi3(1720). The dotted line 
marks the E° threshold, and the dashed line marks the Di3(1895) position. 
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FIG. 11: Differential cross sections for 7 + p — > if + + E°. The number in each panel designates W* (= The solid lines are 
results of the amplitude fits (Eq.[3J discussed in the text. 



the S° angular distributions are shown in Fig. El Com- 
paring the A to the TP shows that in the E° case the 
D wave amplitude plays a more important role, falling 
and rising with a centroid near 1.85 GeV. The P wave 
shows no strong bump in the E°, unlike the A. In this 
case, the G wave coefficient is statistically significant but 
shows little structure. For completeness, we also show 
the magnitude-level fit according to Eq.QJin Fig. El The 
coefficient C\ shows some structure, again due to 5 — P 
or higher-wave interference. The coefficient C2 clearly 
falls and rises, which can be due to P wave activity or 
interferences between S and D waves, for example. 

Figures El an d El show selected differential cross sec- 
tions from this experiment compared to previous data 
and with three published model calculations. The se- 
lected panels show about 1/6 of our data, in increments 
of AW « 80 MeV to show the trends in the cross sections 
and the calculations; the exact W values were chosen to 
emphasize available comparison data. 

The results for the angular distributions of photopro- 
duction of S° are shown in Figs. 1161 and IT7I Again, the 
panels are selected to increase in steps of about 80 MeV 
in W, also to allow comparison to previous data. 



B. W Dependence 

Resonance structure in the s-channel should appear 
most clearly in the W dependence of the cross sections. 
In Fig. El we show the K + A cross section at selected 
angles. The corresponding information for the K + TP 
channel is shown in Fig. 1191 We discuss these results in 
the next section. 

The full set of numerical results from this experiment 
are available from various archival sources, including a 
Ph.D. thesis 0, the CLAS online database [23, or pri- 
vate communication |28| . 



VI. DISCUSSION 
A. Comparison to Previous Data 

FiguresElElEl and !17l show a sample of the differen- 
tial cross section for A and S° hyperon photoproduction 
as a function of angle for a set of W values. For com- 
parison, we can examine the previous large-acceptance 
experiment from SAPHIR at Bonn 0, . There is also a 
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FIG. 14: (Color online) Angular distributions for 7 + p — > + A for selected bins of total energy W . The present CLAS 
results (blue circles) are shown with statistical and yield-fit uncertainties. Data from SAPHIR [Isj (open red stars) and from 
older experiments |40| (black squares) are also shown. The curves are for effective Lagrangian calculations computed by Kaon- 
MAID 5] (solid red) and Ireland et al. Hal (dashed black), and a Regge- model calculation of Guidal et al. |2Cll2l| (dot-dashed 
blue). 
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FIG. 15: (Color online) Angular distributions for 7 + p — > K + + A for selected bins of total energy W. The present CLAS 
results (blue circles) are shown with statistical and yield-fit uncertainties. Data from SAPHIR [lS| (open stars) and from 
older experiments |40| (black squares) are also shown. The curves are for effective Lagrangian calculations computed by Kaon- 
MAID 5] (solid red) and Ireland et al. Hal (dashed black), and a Regge- model calculation of Guidal et al. |2Cll2l| (dot-dashed 
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FIG. 16: (Color online) Angular distributions for 7 + p — > + E° for selected bins of total energy W. The present CLAS 
results (blue circles) are shown with statistical and yield- fit uncertainties. Data from SAPHIR (open stars |1S|] and triangles 
and from older experiments .40.) (black squares) are also shown. The curves are for effective Lagrangian calculations computed 
by Kaon-MAID || (solid red) and Janssen el al. flG^I (dashed black), and a Regge- model calculation of Guidal et al. |20l l2l| 
(dot-dashed blue). 
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FIG. 17: (Color online) Angular distributions for 7 + p — > K + + E° for selected bins of total energy W. The present CLAS 
results (blue circles) are shown with statistical and yield- fit uncertainties. Data from SAPHIR (open stars |1S|] and triangles 
and from older experiments |4(J (black squares) are also shown. The curves are for effective Lagrangian calculations computed 
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FIG. 18: (Color online) Energy distributions for 7 + p — > A" 1 " + A for selected c.m. kaon angles. CLAS results (blue circles) are 
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FIG. 19: (Color online) Energy distributions for 7 + p — > K + +E° for selected cm. kaon angles. CLAS results (blue circles) are 
shown with statistical and yield- fit uncertainties. Data from SAPHIR (open stars [Ts|. triangles 0) and older experiments |40| 
(black squares) are also shown. The curves are for effective Lagrangian calculations c omp uted by Kaon-MAID |g| (solid red) 
and Janssen et al. |Tc|l (dashed black), and a Regge- model calculation of Guidal et al. |20t l2l| (dot-dashed blue). 
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set of measurements that was accumulated from the late 
1950's to the early 1970's using small-aperture magnetic 
spectrometers at CalTech [29ll3(ill3lll3^ . Cornell jsEI^IL 
Bonn IH H3, Orsay |p, DESY|2g, and Tokyo |3|. 
These results are compiled, for example, in Ref [40j. 

The agreement with data from SAPHIR is fair or good, 
but there are some discrepancies. The CLAS results are 
generally more precise, having statistical uncertainties 
that are about 1/4 as large, with about twice as many en- 
ergy bins. The SAPHIR experiment had better backward- 
angular coverage at low energies as well as coverage at 
extreme forward angles where CLAS has an acceptance 
hole. The measurements agree within the estimated un- 
certainties at some angles and generally near threshold 
energies, but CLAS measures consistently larger K + A 
cross sections at most kaon angles and for W > 1.75 GeV. 
This is discussed in more detail below in the context of 
the total cross sections, where it appears that there is 
an energy-independent scale factor of about 3/4 in going 
from the CLAS to the SAPHIR K + A results. The data 
for the K + TP channel are generally in better agreement 
overall: the two experiments agree within their stated 
systematic uncertainties. 

We collected the historic (pre-1973) results from differ- 
ent measurements and plotted them together. The error 
bars are taken as the quoted random uncertainties, with 
no consideration of the quoted systematic uncertainties. 
While these early experiments did not span the large W 
and angular range of the recent experiments, they did 
make high-precision measurements at selected kinemat- 
ics. There are 144 K + A points and 57 K + YP points 
that, overall, are in fair agreement with the CLAS re- 
sults. At backward angles the historic data are in very 
good agreement with the present results from CLAS; at 
forward angles the agreement is fair or good. In the mid- 
range of angles, the historic results are lower than our 
results, and more similar to the SAPHIR data. 

The fit coefficients presented in Figs. [§] and ^| are 
in good qualitative agreement with results published 
by SAPHIR, apart from an arbitrary overall change in 
sign. The CLAS results generally have finer binning 
and smaller estimated uncertainties away from threshold. 
However, our vertical scales do not agree with SAPHIR, 
though it is clear their units are incorrect as given, since 
they should be y^tfe. 

Total cross sections, atot, f° r 7 + P —* K + + A and 
7 + p — ► K + + S° can be calculated from the integrated 
angular distributions. There is some danger in the in- 
tegration procedure since (i) it requires some model of 
the reactions which may bias the resulting fit, and (ii) 
in the absence of complete angular coverage there is also 
the problem of extrapolating the fit into the unmeasured 
section of phase space. Our procedure for extracting and 
calculating the total cross sections was based on fitting 
da/dcos(9'^) in two ways: using Eq.0]to fit the magni- 
tude directly, and Eq.|3to fit the partial wave amplitudes. 
In the magnitude fit, one of the coefficients directly gives 
atot and its associated error. In the amplitude fit atot is 



easily computed from the set of fit parameters, but the 
error is difficult to extract since the fit parameters and 
their errors are correlated. We estimated the systematic 
bias in our integrations by taking the standard deviation 
of the two resultant values as an additional uncertainty, 
and this was added in quadrature to the other estimated 
uncertainties. 

The total cross section results are shown in Figs. |201 
and The error bars combine statistical and esti- 

mated systematic uncertainty due to the fitting proce- 
dure. The gaps in the spectra at W = 2.375 and 2.400 
GeV stem from photon tagger failures at those energies. 
For comparison we show two previously published data 
sets from Bonn 0,^3 S3- Also, bubble chamber data 
for the total cross sections came from Erbc et al. (ABB- 
HHM) 42] . Also shown are model curves for two calcula- 
tions, the effective Lagrangian model embodied in Kaon- 
MAID p| and the Reg ge model of Guidal, Laget, and 
Vanderhaeghen [20J, |2l| . The CLAS results for cr tot differ 
from the Bonn results in an unexpected way, namely that 
the Bonn K + A cross section is smaller than the CLAS re- 
sult by a factor of close to 3/4. This is in contrast to the 
K + TP results, where the CLAS and the Bonn results are 
in good agreement: the values of atot agree well within 
their quoted systematic uncertainties. We note that the 
CLAS results for the two hyperons used exactly the same 
photon normalizations, and that the hyperon yield ex- 
tractions for both cases were made together, as discussed 
above. The acceptance calculations for the CLAS results 
used the same software as well, differing only in the input 
events used for the calculations. In short, we have not 
found any reason within the CLAS analysis for one chan- 
nel agreeing well with previous work and the other not. 
Both results are consistent with the ABBHHM data |43 • 

The CLAS atot results for K + A show a prominent peak 
centered near 1.9 GeV. It does not resemble a simple sin- 
gle Lorentzian, reflective of the expectation that several 
resonant structures are present in this mass range. The 
peak near 1.7 GeV is consistent with contributions from 
the P u (1710) and P 13 (1720). In the case of K+Y, , the 
atot curve shows the previously seen strong peak centered 
at 1.88 GeV, and in addition there is a slight shoulder 
at about 2.05 GeV. The location of the strong peak is 
consistent with the mass of several well-established A 
resonances which may contribute to an isospin 3/2 final 
state. 



B. Comparison to Reaction Models 

The model calculations shown in this paper were not 
fitted to the present results. The effective Lagrangian 
calculations, in particular, were fitted to the previous 
data shown in this paper, and have, therefore, at least fair 
agreement with those earlier results. However, since in 
the case of K + A production we have some disagreement 
with the SAPHIR data in the mid-range of angles, we 
cannot expect these calculations to be in quantitative 
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FIG. 20: (Color online) Total cross section for 7 + p — > A" + + A. The data from CLAS (blue circles) are shown with combined 
statistical and fitting uncertainties. Also shown are results from two publications from SAPHIR (red stars (2004) red 
triangles (1998) |8j), and the ABBHHM Collaboration (light blue squares) 1-1 2B . The curves are from a Regge model (dashed 
blue) Kaon-MAID (solid red) @, Kaon-MAID with the Z?i 3 (1895) turned off (dotted red), and Saghai et al. (dot-dashed 

black) ||. 



agreement with us. It is nevertheless interesting to see 
what the more copious CLAS results seem to indicate in 
comparison to a few of these previous models. 

The Regge- model calculation [23, EI shown in the pre- 
ceding figures uses only K and K* exchanges, with no 
s-channel resonances. The model was constructed to fit 
high-energy kaon photoproduction data 0], f° r W be- 
tween 5 and 16 GeV, and may be expected to reproduce 
the average behavior of the cross section in the nucleon 
resonance region. However, extrapolated down to the 
resonance region, the model overpredicts the size of the 
A cross section and underpredicts that of the S°. This is 
evident in all the graphs, but is especially easily seen in 
the total cross sections, Figs.l2TJlandl2*Tl Since it is a pure 
i-channel reaction model, it cannot produce a rise at back 
angles as seen for the A, and illustrates the need for s- 
and u-channel contributions to understand that feature. 

Two hadrodynamic models 0, Hc| based on similar 
effective Lagrangian approaches are also shown. Both 
emphasize the addition of a small set of s-channel reso- 
nances to the non-resonant Born terms, and differ in their 
treatment of hadronic form factors and gauge invariance 
restoration. As both were fitted to the previous data 
from SAPHIR 0], they are expected to be in somewhat 
poorer agreement with our K + A than our with K + TP 
data. 

Both models contain a set of known s-channel N* reso- 



nances: Su(1650), Pii(1710), andPi 3 (1720). The model 
of Mart et al. f|] which is used in the Kaon-MAID cal- 
culations contains an additional Z?i3(1895) resonance in 
its K + A description. In the K + A case, the calculations 
of Ireland et al. [12^ are shown since they represent an 
update of the earlier work of Janssen et al. [l(| . These 
calculations included photon beam asymmetry |l3j and 
electroproduction data points in the dataset used for 
fitting. The curves displayed on Figs. 1141 1151 and IT51 
contain the set of known resonances plus an additional 
Pn(1895) resonance. This combination was found to give 
the best quantitative agreement with the dataset used for 
fitting. The analysis of Ref. ^3 was restricted to a study 
of the K + A channel, so for comparison with t he p resent 
K + YP data, we use slightly older calculations [T(| which 
contain an additional Z?i3(1895). 

The CLAS K + A results, which show a structure that 
varies in width and position with kaon angle, suggests an 
interference phenomenon between several resonant states 
in this mass range, rather than a single, well-separated 
resonance. This should be expected, since several N* res- 
onances with one- and two-star PDG ratings occupy this 
mass range. From Fig. 1181 the best qualitative modeling 
of the structure near 1.9 GeV at backward angles is given 
by Kaon-MAID [5| , but the model seems to diverge from 
the trends of the data at forward angles. The calculation 
of Ref. 0] gives a poor description of the data in the 1.9 



22 



7 + 



k + + r 



• CLAS (05) 

* SAPHIR (04) 
a SAPHIR (98) 

ABBHHM 




1 .6 1.8 2 2.2 2.4 7.6 

Baryon Mass, W (GeV) 



FIG. 21: (Color online) Total cross section for 7 + p — > if + + E°. The data are from CLAS (blue circles) are shown with 
combined statistical and fitting uncertainties, Also shown are results from two pub lications from SAPHIR (red stars (2004) |T^I . 
red tr iang les (1998) [8j), and the ABBHHM Collaboration (light blue squares) |42|. The curves are from a Regge model (dashed 
blue) Hatll and Kaon- MAID (solid red) 0. 



GeV region at backward angles, but at forward angles 
it is similar to the Kaon-MAID calculation. Using the 
model curves as a guide, we see that a fixed position for 
a single isolated resonance near 1.9 GeV is not consistent 
with the small (~ 50 MeV) variation with angle of the 
feature seen in the cross sections. 

In the S° case there is some indication of a structure 
above the large peak at 1.9 GeV between 2.0 and 2.1 
GeV. This shoulder or small bump in the cross section, 
seen in Fig. EH and in the total cross section Fig. is 
not reproduced by either of the hadrodynamic reaction 
models. 



C. Phenomenological i-Scaling 

The forward peaking of the K + A cross section sug- 
gests that there is substantial contribution to the reac- 
tion mechanism by ^-channel exchange, even in the nu- 
cleon resonance region. To test this idea, the data can be 
cast into the form of da/ dt vs. —t, where t is the Man- 
delstam invariant that gives the 4-momentum squared of 
the kaonic exchange particle(s). The conversion of the 
cross section was done using 

da da 1 . . 

H ~ dcos C K ™- X 2k^ ^ ' 



where k is the center of mass momentum of the incoming 
photon and q is the center of mass momentum of the 
produced kaon. In the simplest Regge picture involving 
the exchange of a single trajectory, the cross section can 
be written as |44| 

j / \ 2a(t)-2 

where D(t) is a function of t only, s is a baryonic scale 
factor taken to be 1 GeV 2 , and a(t) is the Regge trajec- 
tory itself that describes how the angular momentum of 
the exchange varies with t. At our kinematics for small 
\t\ we find a(t) « 0, so the leading behavior of the cross 
section is that it approximately scales with s 2 . 

The cross section da/dt for K + A production is plotted 
in Fig. To obtain sufficient statistical precision, bands 
of width 200 MeV were combined as weighted averages 
(amounting to groups of 8 of our actual bins). The lowest 
band, for E 1 = 1.05 ± 0.10 GeV, starts 40 MeV above 
the reaction threshold. We observe in the figure how the 
cross section values fall on smoothly-varying contours as 
a function of — t. There is an inflected fall-off from the 
minimum —t that is similar for all photon energy bands, 
but as \ —t\ increases the fall-off flattens and then becomes 
a rise. Fig. 1231 shows the cross sections scaled by s 2 , and 
it is seen that there is a clear indication of a locus D(t) 
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describing the data over a range of — t. We interpret the 
departures from this locus as the onset of the s- and u- 
channel contributions to the reaction mechanism. At a 
given value of —t the residual spread of the points can 
be used to determine a(t) for this reaction; this work is 
in progress and will be published separately. 

Examination of Figs. [221 and [2B1 shows a progressive 
flattening of the slope in the cross section as | — t\ — > 0. 
This same "plateau" phenomenon was seen in data from 
SLAC 43] taken at E 7 = 5,8, 11, 16 GeV, that is, well 
above the energies of the present results. In the model of 
Guidal, Laget, and Vanderhaeghen [2(j, this effect arose 
from the interplay of degenerate K and K * Regge trajec- 
tories and the requirements imposed by gauge invariancc 
in the model. The intercepts of these trajectories are at 
a(0) = —0.20 and +0.25, respectively, so their average 
is indeed at about 0, leading to the observed s 2 scaling. 
We note that this plateau effect persists well into the nu- 
cleon resonance region, which suggests the importance of 
K and K* exchange throughout this kinematic region. 

The cross section da/dt for the E° channel is shown in 
Fig- El I n this case, the data do not fall in monotonically 
shifting contours as E 1 increases, as was the case for the 
A in Fig. [22 Instead, a more nucleon-resonance domi- 
nated picture is suggested by the crossing of the bands 
of data points. This is emphasized again in Fig. 1251 that 
shows the s 2 scaled cross sections, which in this case do 
not form a tight band of points. There is no consistent 
trend toward a flattening of the slope, as was the case 
in K + A production; in the previously cited theory [2£j 
this is because in K + T,° production the K plays little 
role compared to K* since gKT,N < 9k an- Further- 
more, the large "resonant" rise in the S° cross section 
near W = 1.90 GeV is serving to cover up any simple 
t-channel behavior for this hyperon. 

At high enough energies, it is expected, however, that 
the E° cross section should also behave as expected by 
t-channel dominance. In Fig. [2EI we show the subset of 
the data from the previous figure for E 1 > 2.39 GeV, 
where the scaling by s 2 does seem to work. We note that 
this is well above the large "A" peak in the total cross 
section, Fi g. |21 J and spans the range where the Regge 
calculation [2(J, * s successful in explaining these data. 

VII. CONCLUSIONS 

In summary, we present results from an experimental 
investigation of A and S° hyperon photoproduction from 



the proton in the energy range where nucleon resonance 
physics should dominate. We provide the to-date largest 
body of data for these reactions in coverage over energy 
and meson angle. Our K + A cross section results reveal 
an interesting lU-dependence: double-peaked at forward 
and backward angles, but not at central angles. We see 
that the structure near W = 1.9 GeV shifts in position 
and shape from forward to backward angles. This finding 
cannot be explained by a i-channel Regge-based model 
or by the addition of a single new resonance in the s or u 
channel. The S° results confirm a single large maximum 
in the cross section near 1.9 GeV, with weak indications 
of more structure between 2.0 and 2.1 GeV. The results 
are in fair or good agreement with several older experi- 
ments. For the K + A case we see that a phenomenological 
scaling of the i-dependence of the cross section by s 2 is 
quite successful in describing the full range of forward- 
angle data, and that this scaling does not work as well 
for the K + Y,° data. Our results show that hyperon pho- 
toproduction can reveal resonance structure previously 
"hidden" from view, thereby improving our understand- 
ing of nucleonic excitations in the higher mass region 
where data are sparse. Comprehensive partial wave anal- 
ysis and amplitude modeling for these results can there- 
fore be hoped to firmly establish the mass and possibly 
the quantum numbers of these states. 
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FIG. 22: (Color online) The entire 7 + p — + + A data set shown as da/dt vs. — t, for ten bands of photon energy with 
AEy — 0.20 GeV. No scaling was applied. 
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FIG. 25: (Color online) The entire 7 +p — > K + + E° data set shown as da/dt vs. —t, for ten bands of photon energy with 
AE~, = 0.20 GeV. The cross sections were scaled by s 2 = W 4 , showing a less well-defined band of points than in the K + A case. 
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